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Abstract 

In this paper we define a causal Lorentz covariant noncommutative (NC) classical Electrody- 
namics. We obtain an explicit realization of the NC theory by solving perturbatively the Seiberg- 
Witten map. The action is polynomial in the field strenght F, allowing to preserve both causality 
and Lorentz covariance. The general structure of the Lagrangian is studied, to all orders in the 
perturbative expansion in the NC parameter 9. We show that monochromatic plane waves are 
solutions of the equations of motion to all orders. An iterative method has been developed to solve 
the equations of motion and has been applied to the study of the corrections to the superposition 
law and to the Coulomb law. 
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1 Introduction 



Suggestions on the possibility that Nature could allow for noncommuting spatial coordinates, came 
both from the past Q and more recently in the realm of superstring theory studying low energy 
excitations of D-branes in a magnetic field || . 

This has stimulated investigations on the noncommutative (NC) versions of gauge field theories 
and on the behaviour of their quantized counterparts. Among these, Maxwell theory is perhaps 
the easiest example and one where a possible experimental test of this hypothesis could be realiz- 
able. Nevertheless, two main problems arise when one tries to implement Electromagnetism in a 
noncommutative geometry: the loss of causality due to the appearence of derivative couplings in 
the Lagrangian and, more fundamentally, the violation of Lorentz invariance exhibited by plane 
wave solutions [H . These problems have been discussed with a different approach in the framework 
of NC QED ||f|. 

In this paper we show that both these problems may be avoided if one allows a nonzero 
" electrical" component into the tensor 9 of the noncommutation relations so including time as a 
NC coordinate. After application of the Seiberg-Witten map || the theory is perturbative in 9 and 
classical plane waves turn out to be exact solutions. They no longer obey a superposition principle. 
Finally a sort of Electric-Magnetic duality comprehending 9 and reminiscent of the known one in 
commutative Maxwell theory, appears between the fields in the equations of motion reinforcing the 
interpretation of 9 as a sort of background primordial Electromagnetic field. 

In Section 2 we fix notations and conventions, recall the definition of the S-W map j2| and 
show the explicit solution to second order in 9. 

In Section 3 we prove that the Lagrangian of the theory is polynomial to all orders in the 
perturbative parameter so that causality is preserved. The equations of motion are derived in 
Section 4 where evidence is also given of the mentioned duality. 

In Section 5 a general iterative method of solving the equations of motion is outlined. 
After proving that plane waves are solutions, the method is applied to the problem of plane wave 
superposition and to derive corrections to the Coulomb law. 

The paper ends with some comments on the results found and on possible experimental settings 
aimed to directly measure noncommutativity. 

2 The S-W map and second order expantions 

In the following, a hat over a classical symbol will indicate the same quantity in its NC version. In 
this fashion, coordinates of flat noncommutative Minkowsky spacetime will be assumed to satisfy: 

[x^,x"}*=i9^ (2.1) 
where 9^ v is a real skew tensor whose components are set as follows: 

QQi 

0ij = e ijkp k ( 2 - 2 ) 

Note that we do not impose e l = 0. This means that time does not commute with spatial coordi- 
nates and 9 is a constant tensor field. Besides 9 we consider the usual Electromagnetic field whose 
NC action is given by: 

S = ~\ j d*xF> lv *Fr, = -~ Jd 4 xF^F^ (2.3) 

The corresponding Lagrangian and field strenght are given by: 

t = -\P^F^ (2.4) 
F„ v = d»A v - cU M - i (2.5) 
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Here the star-product (*) between NC quantities is defined as usual: 

(/ * g)(x) := e^" d ^f{x)g{xX =x , (2.6) 
Also, the following conventions will be used for Electromagnetic fields: 

E l = F oi 

B k = he ijk FV ( 2 - 7 ) 

Now, according to Seiberg and Witten ^ , every NC gauge theory A M has a perturbative description 
in terms of the non commuting parameter 9 and another commutative theory possessing the 
same degrees of freedom as the NC one. The relation between them is established by means of the 
Seiberg- Witten map: 

/ & = -|{ia,^i M + F^}*-(a^/3) (28) 

{ A^g=o = A^ 

Solving the above equations means determine each piece of the perturbative expansions: 

K = Ap+AM + AW + ••• (2.9) 
h v = F MV + F$ +£$+.-. (2.10) 

relating at every order in 9 the NC quantities with their respective classical counterparts. As is 
well known ||, one obtains to first order: 

4 13 - -h^MfyAn+Fpn) (211) 

F$ = 9 jS (F^Fyg — AydsF^) 

Considering the second order corrections, we assume Affl = \ 9 a/s 9 jS n fia p 1 s and substitute the 
whole expansion of A^ into (2.£). We realize that differentiating and then evaluating at 9 = 0, 
we end with a recursive relation between second order and first order corrections and their 9- 
derivatives. This leads to computation of the term n. After careful rearrangements, the expression 
for the second order correction to A^ is: 

A (2) = 1 gafiffys ^ dgAadfiAii + A^FscFfr + A a A 7 d s F Pli + i d^(A a A 7 d s A )^ (2.12) 

Similarly, via the relation F$ = d^Al 2) + 9i 5 d 1 A^ ) d s A v -{jx^v) one also computes: 

p(2) = eaf 3 e js FmFsaFpy _ grrS Aj d s F$ - ~ 6^6^ A 1 (d a A s + A a d s ) d p F^ (2.13) 

3 The general structure 

We discuss some properties valid to all orders in 9 of the perturbative action obtained by means 
of the S-W map. 

Proposition 3.1 The Lagrangian L corresponding to the action ( |£.^ ) via the S-W map is a poly- 
nomial in F only (that is: it does not contain derivatives of F); furthermore the terms £A n ) of 
order n in 9 form a homogeneous polynomial of degree n + 2 in F. 

Proof : 

From the S-W equation (2.S) we have: 



SF^ _ 1 

86<*P 8 ' 



■„ [AcdpA,, + - (a« «-» v) + 1 {A a> dpAf, + fyp} , A t 
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A^{A a ,dpA v +Fp v \] -(a^[3)+S* (3.1) 



Here 5* is supposed to include all the terms arising whenever the derivation acts on the 9 s appearing 
in the * of the * - product; they always give rise to total derivatives in the Lagrangian density and 
so may be neglected. As a consequence, performing an arbitrary number of derivations and then 
putting 8 — shows that commutators of the type present in (|]l]) give vanishing contributions. 

Then all significant contributions arc seen to come from the term i d v ^A a , dpA^ + P^u\ which, 

evaluated at after k derivations, produces an homogeneous polynomial of order k + 1 in A with 
k + 1 derivatives (with respect to spacetime coordinates) equally distributed on each monomial. 
Finally, considering F^" * at order n in 9, by the same argument, one obtains an homogeneous 
polynomial of order n + 2 in A with n + 2 spacetime derivatives comparing in each monomial. 
Now, since the Lagrangian density (obtained from the S-W map) is certainly invariant under the 
usual U(l) gauge transformations, every monomial can be rearranged, modulo integration by parts, 
so as to depend only on F and possibly its derivatives. But being the number of derivatives exactly 
equal to the number of As in every monomial, it follows that derivatives of F cannot appear at 
all. QED 

Corollary 3.1 The equations of motion of the U(l) theory take the form: 

d u F» v = (3.2) 

where F^ u is the sum of homogeneous polynomials of degree n + 1 in F and order n in 9 (i.e. 
written symbolically): 

F = ^9 n F n+l (3.3) 

n 

As we will see this property helps to derive a recursive algorithm for their resolution. 



The main consequence of the structure (3.3) evidenced above, is that the equations of motion 
for the field strenght are of first order. This seems to suggest that the theory is causal even 
though not requiring time commutativity. In the literature, it is suspected that causality does not 
survive Noncommutativity fl. In our model though, after undertaking the SW map, the action 
has been manipulated and integrated by parts to render all terms explicitly gauge invariant. As a 
by-product, all higher order time derivatives have disappeared. In effect, this task is equivalent to 
add boundary terms to the Lagrangian: exactly those capable of giving causal consistency to the 
theory. 

Probably this should be the right procedure to follow generally. Furthermore, the fact that 
preserving causality is no more consequence of imposing zero temporal components in 9, allows 
to require that it can transform like a tensor in respect to the Lorentz group. It descends that 
Lorentz covariance is also preserved. 

4 Equations of motion up to second order 

Let us expand also the NC Lagrangian density fl2.4| ) into pieces of increasing order in 9: 

t = £ + ,C (1) +£ (2) + ••• (4.1) 

The first term here coincides with the classical Maxwell Lagrangian while the other terms are its 
various corrections. More precisely: 

r — _i f 



1(2) = _I{^ 1 )FW+2F^FW} 
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Recall j7J that up to first order in 6, the NC Lagrangian has the following form: 

1 



I - 2 0afjF ^ ) F ' 2 + 20 a/, F> a 'F lux F vP 



or, upon substitution according to our conventions (|2.2|) we have: 

£=i(l + 0- B- e- E)(E 2 — B 2 ) — (/3 • E + e • B)(E • B) 

Next, looking for the second order term, one finds with a little effort: 
1 



(4.2) 



(4.3) 



0Q /9 07 i {^f^ + 2F^F^F^F Sa + F^F^F^F^ 



4 

■ gFp a Fs~fF 2 + -FpjFagF 2 } 



(4.4) 



Here again, after substitution and accurate computation, you get: 



£< 2 > 



l 



(e ■ E — /3 ■ B)(/3 ■ E + e ■ B)(E • B) + - [(e • E — /3 • B) 2 (E 



+ (e • /3)(E 2 - B 2 )(E • B) - (E • B) 2 (e 2 - /3 2 ) ] 
As already remarked, the second variation of £> yields the usual equations of motion: 

d u F^ = 



(4.5) 



(4.6) 



and eq.Q leads us to write F = F + F^ + F^ + ■■■ where F^ = d n F n+1 . 

It is now tempting to regard each piece like this as a correction to the classical field straight 
F due to NC geometry. This is more properly done here than on the expansion ( 2.1 0| ) because 
we are referring to the equations of motion. Furthermore the interesting thing || is that denoting 
the content of the NC field F with an Electric displacement and Magnetic induction (D,H) and 
restating the above expansion as: 

D = E + D (1) +D {2) +■■■ 
H = B + H (1) +H (2) + •■■ 

where the classical fields (E,B) in F are recaptured as their zeroth order corrections (D*-°' ) , H^ - 1 ), 
then the equations of motion take the usual Maxwell form: 



dB 

~dt 

<9D 

~dt 



V x E = 
V-B = 

V x H = 
V-D = 



(4.7) 



(4.8) 



Note that the first two are simply the Bianchi identities; the other two really describe the behaviour 
of NC Electromagnetism in empty space. Working with the first order correction to F which is: 



we obtain for the NC fields the approximated expressions: 

D = (l+/3-B-e -E)E-(/3-E + e-B)B-i(E 2 -B 2 )£-(E-B)/3 



H 



(1 + p ■ B - e • E) B + (p ■ E + e • B) E - i (E 2 - B 2 ) (3 + (E • B) e 



(4.9) 



(4.10) 
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Here the NC ten sor 8 has been assumed to represent a couple of fields (e , P ) in agreement with 



the conventions (2.2). The second order correction to F reads explicitly: 



4 

+ I {e al3 e lS F^F u p F aS + e al3 F^F 5a {9 vS F^ - 0^F V ~<) + 9^9~" J F Sa F Sl F a p} 

- \ {{6F){0^ F%F" + e v(s F^F- lfj - ff^F^F^p) + 6^9^ F aS F Sl F aP ) 

+ \{6F){9 tlv F' 2 + {6F)F» V \ 
8 

+ i (8 ati 9 vS F a 5F 2 + ff^^FasF^F^) (4.11) 

This rather involved formula, when re-expressed in terms of the classical fields gives us the second 
order terms in 9 to be added to the above: 

D (2) = [(e ■ E - (3 • B) 2 - e 2 B 2 + (e ■ B) 2 + (e • /3)(E • B)] E 

+ [((3 ■ E)(e • E - ■ B) - (e • B)(/3 ■ B) + /3 2 (E • B) + 1 (e • /3)(E 2 - B) 2 ] B 

+ [(e- E- p ■ B)E 2 + (/3 • E)(E ■ B) + (p ■ B)B 2 ] e 

+ (e- E- P ■ B)(E • B)/3 + [E- (e x B)]/3 x E (4.12) 
while for the magnetic induction we get: 

H (2) = [£ 2 (E-B)-(e-B)(e-E-/3-B)-(e-E)(/3-E)-i(e-^)(E 2 -B) 2 ]E 

+ [(e • E - p • B) 2 - /3 2 E 2 + (P • E) 2 + (e • /3)(E • B)] B 
- (e ■ E - /3 • B)(E • B) £ + [B • (/3 x E)] e x B 

+ [(£ • E)E 2 + (e • B)(E • B) - (e • E - p ■ B)B 2 ] p (4.13) 

We end this section observing that applying an Electric-magnetic duality directly on the classical 
fields and reversely on the noncommuting parameter in this way: 

E B (414) 



induces, up to second order, an "Electric-magnetic duality" on the NC fields (4.10): 

-D > H , 
H -> D < 4 ' 15 ) 

At present, the meaning of this symmetry is unclear and we suspect it remains true to all orders 
in the perturbative 9 expansion. 

5 Exact solutions and an iterative method 

We seek solutions to the equations of motion: 

d[vF^ p ] = 

(5.1) 

d v F^ = 
where: 

F = F + F {1) +F {2) +■■■ (5.2) 
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with the structure = 9 n F n+2 already evident for example in eq. (JO) and (fhlip 



The most natural thing to suppose is that also a solution should be written as a sum: 

F:=F (0) +F (i) +F (2) + ... (5.3) 

with pieces F^ k ' now understood to be corrections to a solution F^- ' to the classical Maxwell 
Equations i.e. d u F^ fJ/l/ — plus the Bianchi identities. We will briefly state this as dF^ = 0. 
Furthermore, let \ k be the operation of keeping, in a generic expression, all terms up to a given 



order fc in 8, neglecting the others. Then extracting fc-th order from (5.2) terms like this: 



F\ k = Y, F(l) (5-4) 

2 = 

will be present. Accounting for that, hypotesis (|]^) and the structure ( |3.3[ ) we get: 

F\ k =F\ k + e(FF)\ k ^ + --- + 9_^(F_^F)\ (5.5) 

fc fe+i 

Our purpose is to write down a recursive method of solving the noncommutative Maxwell equation 
dF = having a classical solution F(°\ This is realized order by order noting that (dF)\k = dF\ k . 
Then taking first order into the recursive relation (p75J) we have: 



dFli = d(F^ + + 0(FF)\ o ) = (5.6) 
Now, being F^ a classical solution, we are led to solve the equation: 

<9F« = -d(9F^F^) (5.7) 
In exactly the same way, solutions correct up to second order come from: 

of^ = -a [e (f^fW + FWr(°)) + ee f^f^f^A (5.8) 

Generally, obtaining the fc-th term in the expansion ( |5.3| ) always reduces to solving an equation of 
the form: 

d vF W» v = J"[FW,... ,F< k - 1 '>] (5.9) 

where the right member J M only involves all the fc — 1 solutions computed in the previous steps. 
Now, deciding that each two form F^ comes from a potential A^ k ' satisfying the Lorentz gauge 
constraint^ d v A^ k ' v = then Eq.(5.9) becomes: 



- DA 1 - 1 = (5.10) 

This is immediately solved employing the Lienard-Wickert potentials. Then in principle we have 
got an authomatic tool capable of solving the equations of motion in full. 
Let us focus, for example, on the single plane wave solution: 

A^^e lk ' x (5.11) 

with fc M fc M = Cuk^ = in the Lorentz gauge d v A v = 0. We have: 

=i{k^v-K^y k -* (5.12) 

This is a particular case because we will now show that it is an exact solution of eq.(|]l]). 

Lemma 5.1 Given an antisymmetric matrix 0^ , a null vector k a and a family of vectors {C^}iei 
orthogonal to k a , then any combination of n copies of 0^ v , (n+1) vectors of the given family and 
(n+2) copies of k a in which all indices but one are saturated, vanish. 

1 It is easy to show that this can always be done order by order 
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Proof. Try to build a nonvanishing combination. In so doing, you cannot saturate the k vectors 
with the C vectors due to ortogonality. Neither you can saturate two of them with one 9 matrix 
due to its antisimmetry. You are obliged to use one only k vector for each matrix, spending n of 
them. Of the two remaining, one can be chosen as the free index but the other must necessarily be 
saturated with one of the £ vectors or one of the 9 matrices giving a vanishing result. QED. 



Proposition 5.1 Monochromatic plane waves solve the field equations (5.1) to every order in 9. 
Proof. Let us write the general monochromatic plane wave as: 

A M = $ M (tf -x) (5.13) 
with K 2 = and K^<f>'^(K ■ x) = so that 

F^=Kp& v {K>x)-K v & ll {K-x) (5.14) 



Let be the term of order n in 6; then d^F^^ is the sum of terms obtained by contraction 

of n copies of 9^ v , n copies of $' Q , one copy of $>'l and n + 2 copies of K 1 . From the Lemma it 
follows that d^F^" = . QED. 

The previous property of monochromatic plane waves holds for any lagrangian having the assumed 
polynomial structure, independently of the fact that it has been derived from a NC theory using 
the SW map. 

5.1 Plane wave superposition 

While single plane waves turn out to be exact solutions of the field equations this is no longer valid 
even for a simple superposition like this: 

A^.= (,e lk - X + Cy k '- X (5.15) 
corresponding to the classical solution (by linearity): 

F$ ^iW-We" 1 +i{k'^' u -k'„C'y k '- x (5.16) 



To find out its first order correction in the NC framework we must solve (5.7) yielding: 

d v F^ = i { [(fcfc'MCC) - (CC)e(kk')} (ft" - O 

- PCMCfc') - «k') e (kC)] (ft^ + ft" 1 ) 
+ 2 [{kC)e(kk') - {kk') e (kC)]C 

- 2 {(Ck')e(kk') - (kk')e(Ck')} C } e l ^ k+k '^ x (5.17) 

where the following anti-symmetric inner product has been defined: (vw)e '■= w M (9 M „ w v . 
This equation can be solved assuming 

Fj$ = - c^W (5.18) 

and still satisfyi ng an extended Lorentz gauge constraint d v A^ v = 0. 
Infact, rewriting eq. ( 5.17 ) in the form 

-□4 1 ' =iJ lx e i( - k+k ">- x (5.19) 

we realize that J is transverse to k + k! as can be easily proved using the defining relations: 

k 2 = ft' 2 = k ■ C = ft' • C' = (5-20) 

This means that if we put abruptely, 

A W = _!^L_ e i(fc+fe')-x (521) 



" 2k -k 
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this solves eq.(5.17) being also compatible with the extended Lorentz gauge. 



Note that this corrected version of the superposition law could be used to reveal a refraction 
effect suffered by a ray of light in passing from an empty region to one in which a background 
static and uniform magnetic or electric field is present j| [|. The incoming and reflected rays 
propagating in the empty region, should be described by a superposition of waves agreeing 
with the refracted one in the transition region. 

5.2 The Coulomb Law 

All NC theories are characterized by a parameter 9 which defines a natural scale of length. 
From a dimensional analysis, the corrections to the Coulomb law are of order but a 
complete power series expansion in is expected so that if L is its convergence ratio 
(plausibely finite), then non-perturbative contributions should become relevant in the 
region r < y/6L where the perturbative description fails. We can make a sensible study of 
the NC corrections to the Coulomb law, considering the potential generated by a charged 
conducting sphere of radius ro- 
At 0-th order the classical potential is 



A (0) 



■■pdt r>r 
~dt r<r 



(5.22) 



so that {x l is the radial versor} 



r 0i 



-^x l 9(r-r ) 



(5.23) 



At first order: 



Q F (1)„U = _Q G „U 



(5.24) 



with 



G?» ■- -\qf)F» v - - 9^ U (FF) + 9 aP F fia F u P + (9^F au - 9 uf3 F af *)F c 



A direct computation gives for the tensor G: 

„2 



so that 



G 0i 

Qij 



(e • x)x l 



[r - r ) 



1 



9(r - r ) 



a 3 



(5.25) 

(5.26) 
(5.27) 



(5.28) 



with 



J 



(4e-r;fA/3) 9(r - r ) + (3/2 e • f; -1/2 r A 0) 5{r - r ) 



(5.29) 



9 




-100 -50 50 100 



Figure 1: The correction to the Coulomb law due to The figure shows the equipo- 

tential levels in the (x, z) plane, assuming that e is oriented along the z-axis. We have 
chosen e = 1, rg = 1. The length unit is r§. 



We have solved numerically the equation 5.28 and we show in Figure [l] the equipotential 
level of the zeroth component of A^> which is symmetrical under rotations about the 
direction of e. The other components have a similar angular behaviour and their precise 
values depend on the direction and magnitude of j3. 

The corrections in Figure |l| give the modification of the Coulomb law in the case of a 
charged sphere. This interpretation would be well defined with ro ^> y/6, being the 
perturbative solution valid almost everywhere, even inside the conducting sphere. But in 
the specific example considered in Figure [l] the sphere coincides with the excluded region, 
where the perturbative approach fails. This case suggests a different interpretation, as the 
NC correction to the Coulomb potential of a point-charge. In fact in NC theories it is 
intuitive to replace pointlike with extended object, whose typical length is s/6. 
The corrections to the potential violate the Gauss law and the spherical symmetry of the 
classical solution. As a consequence, we observe in the case of a conducting macroscopic 
sphere, that the potential inside the conductor is not constant. This remark suggests a 
way to test NC electrodynamics effects. In Figure we show the relative contribution of 
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50 



100 



Figure 2: The ratio /A^> which shows the effect of the NC corrections to the Coulomb 
law, in the (x,z) plane. We have chosen e oriented along the z-axis, e = 1, e = 1 and 
ro = 1. The length unit is tq. 

the corrections A^' to the classical Coulomb potential A^ . The size of the corrections is 
already relevant (e.g. greater than 10 percent) at a length scale bigger by more than one 
order of magnitude w.r.t. to the one determined by the NC parameter. 

6 Conclusions 

In this work we have formulated an explicit perturbative realization of NC electrodynam- 
ics, which turns out to be causal and Lorentz invariant. The basic steps to obtain this 
result have been: the use of the SW map and a rearrangement of the action aimed to 
render every term explicitly gauge invariant, by use of careful integration by parts. The 
resulting expressions do not contain time derivatives of order higher than two, yielding 
authomatically a causal theory. This latter property is obtained without imposing any 
constraint on the NC parameter 6, which can be chosen in full generality as a Lorentz 
tensor, leading to a Lorentz covariant theory. 

We have studied the general structure of the Lagrangian, to all orders in the perturbative 
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expansion. We have shown that the monochromatic plane wave is solution of the equations 
of motion to first Q and even to all orders. 

We developed an iterative method to solve the equations of motion. In particular we 
applied this method to study the corrections to the superposition law of plane waves and 
to the electrostatic potential of a spherically symmetric charge distribution. The most 
relevant qualitative feature of the NC corrections that we calculated is that they have a 
peculiar signature which makes them, at least in priciple, distinguishable from the classical 
corresponding effects. A possible test of the superposition law could be done by studying 
the reflection and rifraction of light on a magnetic field, using for instance the experimen- 
tal setting described in flO|| . Furthermore, the deviations from the Coulomb law could be 
evidenced by measuring the charge distribution on the surface and the electric field inside 
an empty conducting sphere. 
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